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Abstract. In this paper we consider dividend problem for an insurance 
company whose risk evolves as a spectrally negative Levy process (in the 
absence of dividend payments) when Parisian delay is applied. The ob- 
jective function is given by the cumulative discounted dividends received 
until the moment of ruin when so-called barrier strategy is applied. Ad- 
ditionally we will consider two possibilities of delay. In the first scenario 
ruin happens when the surplus process stays below zero longer than fixed 
amount of time > 0. In the second case there is a time lag d between 
decision of paying dividends and its implementation. 
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1 Introduction 

In risk theory we usually consider classical Cramcr-Lundberg risk process: 



where a: > denotes an initial reserve and Ui, {i — 1, 2, ...) are i.i.d distributed 
claims with the distribution function F. The arrival process is a homogeneous 
Poisson process Nt with intensity A. The premium income is modeled by a 
constant premium density c and the net profit condition is then XEUi/c < 1. 
Lately there has been considered more general setting of a spectrally negative 
Levy process. That is, X = {Xt}t>o is a process with stationary and indepen- 
dent increments having nonpositive jumps. We will assume that process starts 
from Xq = X and later we will use convention P(-|Xo = a:;) = and Pq = P. 

Such process takes into account not only large claims compensated by a steady 
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income at rate c > 0, but also small perturbations coming from the Gaussian 
component and additionally (when v(—oo,0) = oo for the jump measure v of 
X) compensated countable infinite number of small claims arriving over each 
finite time horizon. Working under this class of models, it became apparent 
that, despite of the diversity of possible probabilistic behaviors it allows to ex- 
press all results in a unifying manner via the c-harmonic scale function W^'^'^ (x) 
defined via its Laplace transform. This paper further illustrates this aspect, by 
unveiling the way the scale functions intervenes in a quite complicated control 
problem. 

The classic research of the Scandinavian school had focused on determining 
the "ruin probability" of the process ^ ever becoming negative, under the 
assumption that X has positive profits. Since however in this case the surplus 
has the unrealistic property that it converges to infinity with probability one, 
De Finetti |14) introduced the dividend barrier model, in which all surpluses 
above a given level are transferred (subject to a discount rate) to a beneficiary, 
and raised the question of optimizing this barrier. Formally, we consider the 
risk process controlled by the dividend policy tt given by 

U^^Xt~L^, (2) 

where Xq = x > is an initial reserves and is an increasing, adapted and 
left-continuous process representing the cumulative dividends paid out by the 
company up till time t. The optimization objective function is given by the 
average cumulative discounted dividends received until the moment of ruin: 

v^x)^E, f e-^'dL^, (3) 
Jo 

where a'^ is a ruin time that we specify later depending on the considered 
scenario and g is a discounting rate. 

The objective of beneficiaries of an insurance company is to maximize (x) 
over all admissible strategies tt: 

v„{x) — supv^{x), (4) 
Tren 

where 11 is a set of all admissible strategies, that is strategies tt = > 0} 

such that - < C//L. 

For classical risk process ([IJ an intricate "bands strategy" solution was dis- 
covered by Gerber [TS], [T^, as well as the fact that for exponential claims, this 
reduces to a simple barrier strategy: "pay all you can above a fixed constant 
barrier a" . 

There has been a great deal of work on De Finetti's objective, usually con- 
cerning barrier strategies (see e.g. Schmidli [33] for more detailed overview). 
Gerber and Shin 17 and Jeanblanc and Shiryaev |21| consider the optimal 
dividend problem in a Brownian setting. Irback |20| and Zhou |32| study the 
constant barrier under the Cramer-Lundberg model ([T|) . Hallin [19] formulated 
time dependent integro-differential equations describing the payoff associated to 
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a 2n bands policy. The optimality of the "bands strategy" was recently estab- 
lished by Albrecher and Thonhauser [1] in the presence of fixed interest rates as 
well. For related work considering both excess-of-loss reinsurance and dividend 
distribution policies (e.g. in a diffusion setting), see Asmussen et al. |3] and 
[28j and references included in this papers, and for work including also a utility 
function, see Grandits et al. |18j . 

For Levy risk process considered in this paper without Parisian delay Avram 
et al. [5], Kyprianou and Palmowski [53], Loeffen [55] and Loeffen and Renaud 
|27] found sufficient conditions for which barrier strategy is optimal. In fact 
Avram et al. [7] prove that in this case the bands strategy is optimal. 

In this paper we want to analyze the dividend problem where there is so- 
called Parisian delay either at moment of payments of dividends or at the ruin 
time. The name for this delay comes from Parisian option that price are acti- 
vated or canceled depending on the type of option if the underlying asset stays 
above or below the barrier long enough in a row (see [5] and [T5]). 

Dassios and Wu [15] for classical risk process ([IJ consider a Parisian type de- 
lay between a decision to pay a dividend and its implementation. The decision 
to pay is taken when the surplus reaches the fixed barrier a but it is implemented 
only when the surplus stays above barrier longer than fixed d > 0. The dividend 
is paid at the end of this period. This strategy we will denote by tTq. Similar 
problem for a spectrally negative Levy process of bounded variation was ana- 
lyzed [251 . In this paper we generalize this result into general spectrally negative 
Levy risk process. In this case ruin time is given by: a^" — mi{t > : U^'^ < 0}. 
Since the ruin time is classical we know that band strategy is optimal and we 
also know the necessary conditions under which an optimal strategy is the bar- 
rier strategy. We still believe that this new Parisian strategy (although not 
optimal within all strategies) could be very useful for the insurance companies 
giving possibility of natural delay between decision and its implementation. 

In this paper we also consider Parisian delay at the ruin. We denote this 
strategy by 7r°. That is ruin occurs if process If^ stays below zero for longer 
period than a fixed C > 0. Formally, we define Parisian time of ruin by: 

(7"'° = M{t >0:t- sup{s < t : Uf > 0} > C, Uf < 0}. (5) 

We first analyze the strategy tt" according to which the dividends are paid 
according to classical barrier dividend strategy transferring all surpluses above 
a given level a to dividends. We also prove the verification theorem for this type 
of ruin. In particular we find sufficient condition for the barrier strategy to be 
optimal. 

In fact combination of both scenarios is also available. To simplify analysis 
we decided to skip this possibility. 

We believe that giving possibility of Parisian delay could describe better 
many situations of insurance company. For example it can checked if if indeed 
company's reserves increase and we can pay dividends (in the first scenario) or 
it can given possibility for the insurance company to get solvency (in the second 
scenario). 
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The paper is organized as follows. In Section [3] we introduce basic notions 
and notations. In Section [3] we find discounted cumulative dividends payments 
until Parisian ruin time. In Section we prove the verification theorem and 
find necessary conditions for the barrier strategy to be optimal. In Section [5] 
we analyze the case when there is a time lag between decision to pay dividends 
and its implementation. 



2 Preliminaries 

We first review some fluctuation theory of spectrally negative Levy processes 
and refer the reader for more background to Kyprianou [53], Sato [21] and 
Bertoin [8] and references therein. 

In this paper we consider a spectrally negative Levy process X — {Xt\t>o, 
that is a Levy process with the Levy measure v satisfying z/(0, oo) = (for sim- 
plicity we exclude the case of a compound Poisson process with negative jumps). 
Since jumps of a spectrally negative Levy process X are all non-positive, moment 
generating function E[e^'''"*] exists for all > and is given by E[e^"^'] = e*'''*^^-' 
for some function il){9) that is well defined at least on the positive half-axes where 
it is strictly convex with the property that lime_j.oo V'(^) = +oo. Moreover, ip 
is strictly increasing on [$(0),cx)), where $(0) is the largest root of '0(^) — 0. 
We shall denote the right-inverse function of by $ : [0, cjo) — >■ [$(0),(X)). We 
will consider also the dual process Xt ~ ~Xt which is a spectrally positive Levy 
process with the jump measure d (0, y) = v {—y, 0). Characteristics of X will be 
indicated by using a hat over the existing notation for characteristics of X. 

For any 9 for which V'(^) = logE[exp0Xi] is finite we denote by an 
exponential tilting of measure P with Radon-Nikodym derivative with respect 
to P given by 

^ ^exp{eXt-mt), (6) 

where J^t is a right-continuous natural filtration oi X . Under the measure P^ the 
process X is still a spectrally negative Levy process with characteristic function 
ipe given by 

Ms)^Hs + 0)~^P{9). (7) 

Throughout the paper we assume that the following (regularity) condition 
is satisfied: 

.0 

(T > or / xv{dx) = oo or ^{dx) << dx, (8) 



where a a Gaussian coefficient of X. 
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2.1 Scale functions 



For p > 0, there exists a function W^^^ : [0,oo) — >■ [0,oo), called the p-scale 
function, that is continuous and increasing with Laplace transform 



-Ox 



e > <s>{p). 



(9) 



The domain of (p' is extended to the entire real axis by setting W'^p^ (y) = for 
y <0. We denote W'-^^x) = W{x). For later use we mention some properties of 
the function W^^^ that have been obtained in literature. On (0, oo) the function 
y H-)- W^P^{y) is right- and left-differentiable and under the condition (|S]), it holds 
that y I—)- W^P^{y) is continuously differentiable for y > 0. Moreover, if cr > 
it holds that W^p'> £ C°°(0,oo) with W^p'>'{0) = if X has unbounded 

variation with cr = 0, it holds that W'-p^'{0) = oo (see [25] Lemma 4]). 

The function plays a key role in the solution of the two-sided exit 

problem as shown by the following classical identity. Letting t+ , t~ be the 
entrance times of X into [a, oo) and (~oo, —a) respectively, 

T+ = M{t >0:Xt>a}, t- ^ mi{t >Q:Xt<-a} 

it holds for z £ [0, a] that 

-»\ro- > T+] - W^P\z)/W^P\a). (10) 

Closely related to W'^p^ is function Z^p^ given by 

Z(f)(z) = l + gW^^'^z), (11) 

where w''^\z) — W'^P\y)dy is the anti-derivative of W^pK Moreover, the 
scale functions appear also in so-called two-sided downward exit problem: 

M/(p)(z) 



^ ' n 



Z^P\z) - Z'^P\a)- 



(12) 



M^(p)(a)' 

and in one-sided downward exit problem that for any /3 with ■0(/3) < oo> V ^ 
t/'(/3) V and a; > gives: 



where W^"'' and Z^"' are scale functions with respect to the measure P^, u = 
p — ?A(/3) and ^/^(u) is understood in the limiting sense if u = 0. In fact for 
each z G R, W^'p^ (z) is analytically extendable, as a function in p, to the whole 
complex plane; and hence the same is true of Z^p\z). In which case arguing 
again by analytic extension one may weaken the requirement that p > tpil^) V 
to simply p > 0. 

The 'tilted' scale functions can be linked to non-tilted scale functions via the 
(z) = W^P'){z) from O Remark 4]. This relation implies that 



-PT-+I3X _ 
e 



An) 



Tg < oo 



.0^ 



z'-\z) 



$(w) 



(13) 



relation e^^w'"^^ 



Z^^\z) = 1 + u 



-PyW^P\y)dy. 
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2.2 Parisian ruin 

One of most important characteristics in risk theory is a ruin probability defined 
by Vxi^Q < oo) for Tq" = mi{t > . Xf < 0}. Czarna and Palmowski [5] 
extended this notion to so-called Parisian ruin probability, that occurs if the 
process X stays below zero for period longer than a fixed C > (see also [TUl H] 
for the result concerning classical risk process). Let 

= inf {t > ; t - sup{s < t : Xs > 0} > (, Xt > 0} 

and Parisian ruin probability we define as: 

P(r'^ < oo\Xq = x)= Px(r'^ < oo). 

The following result summarize (9, Theorems 1 and 2]. 

Theorem 1 Parisian ruin probability equals: 

P,(t^ =«)) =P,(to- =cx))P(r^ <cx)) (14) 
+ (1 - P(r^ < oo)) (^1 - ^ P(r+ > CWAto < ^, -X^- E dz)^ 

and 

P,(to- =oo) = V'(0+)W^(x), (15) 

nOO POO 

/ e-^" ds / P(r+ > s)Fx{t^ < oo, -X^- E dz) (16) 
Moreover, 

(i) If X is a process of bounded variation, then 

/r > Onro < ^, -X^- E dz) 

P(r^ < oo) 
where 



+ P(r+ > C)P(to- < oo, -X^- E dz) ' 



/•OO 



/ P(r+ > s)F{t^ < oo, -X^- E dz) 
Jo " 

j- r(l-e-*W^)p(z,oo)dz. (18) 



(ii) If X is a process of unbounded variation, then 



V^<oo)=limMM)_^, (19) 
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where 

e-^e-^^qis t) dt ds = (20) 



JO 



m w = lim 7^=^ (21) 

40 n(e) 



/or normalizing function n. 



3 Parisian delay at ruin 

In section we will consider Parisian ruin time ([SJ and dividends paid according 
to barrier strategy that correspond to reducing the risk process U to the level 
a if a; > a, by paying out the amount x — a, and subsequently paying out the 
minimal amount of dividends to keep the risk process below the level a. It is 
well known (see [6]) that for < x < a the corresponding controlled risk process 
U"^ under P^^ is equal in law to the process {a — Yt : t > 0} under P^; for 

Yt = aVXt-Xt (22) 

being Levy process X reflected at its past supremum: 

Xt = sup Xs, 

0<s<t 

where we use notations y V = maxjy, 0}. In this case for all a; > 0, 

Vaix) := v^^ ix) = ^ e-«*dir j , 

and ^ aV Xt — a. 
Note that for x < a, 

Va (x) =. [e-^r^ ^ ^+ < ^C] „^ (a) (23) 

and 

Vaix) = X — a + Va{a) for x > a. (24) 

Assume that X oo a.s. Then by Markov property and fact that X jumps 
only downwards we derive: 

F^ir^ = cx)) = P,(t+ < Ti)Va{T'^ = oo). (25) 

Hence 



\(tC = (X3)- 
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Using change of measure (0) with 9 ~ ^{q), Optional Stopping Theorem and fact 
that on P*('3) process X tends to infinity a.s. (since V'$(g)(0+) — V''(*&('?)+) > 
0), we have for x < a, 



y(9)(a)' 



(26) 



where 



F(«)(x) = e*(«)"P*(')(T^ = oo). (27) 

The probabihty P*^'-'(t^ = oo) hence also function y'^^ could be found using 
Theorem [1] 

It follows from Theorem [T] that under the condition ([8|) function F^'^ (y) 
(similarly like W^'^\y)) is continuously differentiable for y G K. 
Moreover, for n e N, by (IMl) . 



and 



?^a(a) > Ea 

= 

t'a(a) < Ea 



' a+l/n 



' a+l/n 



I a + - 

n 



Va{a.) + - 
n 



' a+l/n 



'Ua(a) 



-E„ 



= Ea 



" a+l/n 



' a+l/?i 



Wa(a) 



1 - Ea 

nq 



' a+l/n 



since LJ^ = Xt — a under can increase only by e up to time 'T'^j^iln Last 
increment in above equation is o(l/n) since by strictly positive drift X is regular 
for (0, oo). Hence, 



F(«)(a) 



1/(9) 



lY Va{a) + - + o - 



and then 



y(9)(c 



y('?)'(a)' 

Thus from (|23)) . (|24p and (p6)) it follows that is continuously differentiable 
for all a; S R and 



V('!)'(o) ' 



X < a. 



(28) 



^ °' + y(9)'(a) ' 2: > a. 
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In particular, 

v» = 1. (29) 
Hence we get the following theorem. 

Theorem 2 The value function corresponding to the barrier strategy tt" is given 
by i28\). The optimal barrier a* is given by: 

a* = inf{a > : V'^''^'{a) < V'^^^' (y) for all y > 0}. (30) 

In particular, if V'^'^'^ G C^(K) and there exists unique solution of equation: 

T/(9)"(a*) = 0, (31) 

then a* is optimal barrier. 

Remark 1 Note that G C^{M.) if W'-'^^ G C^{R). This is the case if e.g. the 
Gaussian component is present. 

4 Verification Theorem 

To prove the optimality of a particular strategy tt across all admissible strategies 
n for the dividend problem (|4]), where the ruin time cr^ is given by the Parisian 
ruin ([5|), we are led, by standard Markovian arguments, to consider the following 
variational inequalities: 

r/(a;)-g/(x) < 0, if x e R, (32) 
f'{x) > 1, if xeR, (33) 

for functions / : M -> R in the domain of the extended generator F of the process 
X which acts on C^{0, oo) functions / as 

r/(a:) = ^f"{x) +p,f'{x) + f [fix + y)- fix) + /'(x)2/l{|,|<ij] i^idy), 

(34) 

where v is the Levy measure of X and denotes the Gaussian coefficient and 
Pq = c — j'^-^yi'idy) if the jump-part has bounded variation; see Sato [311 
Ch. 6, Thm. 31.5]. In particular, if E\X\ < oo and X has unbounded variation 
[bounded variation], a function / that is [C^] on [0, oo) and that is ultimately 
linear lies in the domain of the extended generator. 

Theorem 3 Let C £ (0,oo] and suppose f is continuous and piecewise on 
(— oo,C) if X has bounded variation and that f is and piecewise on 
(— oo,C) if X has unbounded variation. Suppose that f satisfies I13S\) and LSS]) . 
Then f > sup^gjj^^ v'^ for defined in with Parisian ruin time @), where 
Tl<c is a set of all bounded strategies by C . In particular, if C — oo, then 
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Proof We will follow classical arguments. Let tt G TL<c be any admissible 
policy and denote hy L = and U = If^ the corresponding cumulative div- 
idend process and risk process, respectively. By Sato [211 Ch. 6, Thm. 31.5] 
function g{t,x, z) = e~'^*/(2;)l{2<i;} is in a domain of extended generator of the 
three-dimensional Markov process {t, U^, ), with c;^ = t-suj){s <t:Ut>0}. 
Note that finite number of discontinuities in / and hence also single discontinu- 
ity in l{2<f} are allowed here. Hence we are also allowed to apply Ito's lemma 
(e.g. [201 Thm. 32]) if X is of unbounded variation and the change of variable 
formula (e.g. [501 Thm. 31]) if X is of bounded variation: 

e-«*/((70lK-<c} - /(t^o) = Jfit)- f e-^^f'{Us-)dL'^, 

+ f e-''%rf-qf){U,_)ds + Mt, (35) 
Jo 

where Mt is a local martingale with Mq — 0, i"^ is the pathwise continuous 
parts of L and for a function g the process Jg is given by 

Mt) = E e"'' ^9(^^ + - •9(^^)1 l{B.#o}, (36) 

s<t 

where As — Us- + AXs with Axs = Xs — Xs- and Bs ~ —ALs denotes 
the jump of —L at time s. Let T„ be a localizing sequence of M . Applying 
Optional Stopping Theorem to the stopping times T^. = Tk A a'^ and using 
Fatou's Theorem we derive: 

fix) > E,e-''^"/(C/Tjl{,^, <c) - Ji{K)+^. f" e-'^^nUs-)dL's 

Jo 

-E, r" e-i^{rf-qf){Us-)ds. 
Jo 

Invoking the variational inequalities f'{x) > 1 (hence f{As + Bg) — /(^s) < 
-ALs if As > 0) and Tf{x) ~ qf{x) < we have: 

fix) > E,e-«^"/([/Tjl{,.,<^j+E, /" "e-«^di. 











> E, 


'e-'^'fiU.^);a-<K 


+ E, 


/ e-''^dL,;a^ <r; 








Jo 



Letting n — >■ cxi in conjunction with the monotone convergence theorem and 
using fact that l^^u^^^ = complete the proof. □ 
Using verification theorem we find necessary conditions under which the 
optimal strategy takes the form of a barrier strategy. 

Theorem 4 Assume that a > or that X has bounded variation or, otherwise, 
suppose that G C^(0, oo). If q > 0, then a* < oo and the following hold true: 
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(i) 7r° is the optimal strategy in the set IV<a* of all hounded strategies by a* 
and Va' = sup^gn<„. '"'^ ■ 

(a) If {Tva- — qVa-'){x) < for X > a*, the value function and optimal 
strategy of ^ is given by u» = Va" , where the ruin time cr^ is given by the 
Parisian moment of ruin 

The proof of Theorem 0] is based on the verification Theorem |3] and the 
fohowing lemma. 

Lemma 1 (i) We have a* < oo. 

(a) It holds that (Fua* — QVa*){x) = for x < a* . 
(Hi) For X < a* 

<.(x)>l. 

Proof Point (i) follows from the fact that V^^'^-''(j/) is continuous and in- 
creasing starting at some point. Indeed, note that by [22j we have V'-'^'^iy) = 
e'i'(i)yVy{T< = oo) > e*('?)yp^(T-- = oo) = ^t^W^^^?/) and VF(«)'(?/) tends to 
oo as y — > oo. The proof of (ii) follows from (1231) and the martingale property 
of 



I 'a* ^ ' 



= E 



where x < a* . Point (iii) is a consequence of (I29p and definition of a* given in 
(1301). □ 
Moreover, we can give other necessary condition for the barrier strategy to 
be optimal. 

Corollary 1 Suppose that 

^('^'(a) < for all a* <a<b. 

Then the barrier strategy at a* is an optimal strategy. 

Proof Using Theorems [5] and 2] the proof is the same as the proof of }26[ Theorem 
2]. □ 

Corollary 2 Suppose that, for x > 0, v'{x) is monotone decreasing, then tt" 
is an optimal strategy of Q). 

Proof By ((27)) the proof is similar like the proof of [26l Theorem 3] . In fact it 
suffices to prove that V$(q) has completely monotone derivative. This fact follows 

from Theorem [1] since ^P*^''^(to" = oo) and |^P*^''^(to" < oo, -X^- e dz) are 
completely monotone. Indeed, it is known that 

^P*(«)(ro- =(^) =^*(,)(0,0)C/4(,)(0,x), (37) 

where C/$(g) is a renewal function of the descending ladder height process Ht 
under P*^'?) and K<i,(^) (a, /3) is a Laplace exponent of bivariate descending ladder 
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height process {Lt^.Ht) under P*(9) with K<i,(g) (0, 0) = ti^'{^{q)) > 0. From the 
proof of dSl Theorem 3] it follows that U^f^^^{0,x) hence also ■§^^x'"''\t^ — oo) 
is completely monotone. Moreover, by [211 (7.15), p. 195] we have 

px poo 

= (0, 0) / U^(q){x-dy) e-*(«)("+'^)9(dz + w) di;, 
Jq Jo 

and hence ^^x^'^\tq < (X),—X^- e dz) is also completely monotone. □ 



5 Parisian delay at the moment of dividend pay- 
ments 



In this section we analyze the case when we pay dividends only when surplus 
process stay above barrier a longer than a time lag d > 0. The dividends 
are paid at the end of that period and they are paid until regular ruin time 
0-'^" = inf{t > : Ut'' < 0}. Then by (HU]) for x S [0, a], 

= SSM^W; (38) 



and by Markov property for x > a, 

v{x) = e-'^'E^-a [{Xd + t'(a)), To" > d] 

+v{a) [ Ea-y \e"i^^,T+ <T- 

+v{a)E.,^a [e-9^"^X^- = 0,To < d 



iy(«)(a) 

e-^^o,-X^- edy.T^ <d 
(39) 



where E,. 



; ~a ^ "o 



is given in (|T()|) . 

Double Laplace transform of E^ e^'^^" , —X^- E dy,TQ < s for z > by 
13)) equals 



e-"'e-^yE, 



"'[0,00) 



^0 , -X - e dy, To < s 



ds 



= -E, 
a 



e ^0 , Tq < 00 



a V *I'(^g 



■W^^^(^)) :=i?,(/3,^), 



(40) 



where Uq = a + q — tpiP)- Moreover, 



[Xs, T(, > s\ ds ^ — i z - E 



X^-e ,To < 00 



(41) 
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Further, the value v{a) is determined by ([5^ if X has no Gaussian compo- 
nent (ct = 0) or by the smooth paste condition; 



v'{a—) ~ v'{a+), 



(42) 



otherwise. 



Lemma 5.1 If a > then |4^[ ) holds. 
Proof For n G N, 

w(a) = v{a~l/n)Ka 
+v{a + l/n)Ea 



-1/n T ^ + 

' a-l/ri ^ a+l/n 



a+l/n a— l/ri 



1- 



(43) 



where the last term is bounded above by ^{j'^j.^^ > d). Moreover, by ([10 



14^(9) (1/n) 
14^(9) (2/n) 



and by ((T2 

E, 



Z(«)(l/n)-Z(«)(2/n) 



W^'i'>{l/n) 
W(i'){2/n)' 



Multiplying both sides of p3|) by two, subtracting v(a — 1/n) + v{a), dividing 
by 1/n gives: 



v{a) — v{a ~ 1/n) v{a + 1/n) — v{a) 



1/n 



1/n 

v{a + 1/n) — v{a — 1/n) 



1/n 



2VF(9)(l/n) 
M^(«)(2/n) ^ 



(44) 
(45) 



+w(a - l/n)2q 



i:"'w^^Hy)dy~^Smi:^''wi''Hy)dy 



W(iH2/n) Jo 

1/n 



(46) 



where we use ([TT|) . Now, since 

W(9)(0) = and W(9)'(0) = 4, we have: 
VF(«)(1/") M^(''''(0) 1 



lim 



lF(«)(2/n) 2VF(9)'(0) 2' 



Hence increment dH]) converges to u(a)2g(VF(9)(0) - jW'''^\0)) = as n — >■ oo. 
Moreover, 



hm n 

n— >oo 



2Ty('?)(l/n) ^ ^ 
W('!){2/n) " 



1 2/n W^(«)(2/n) - 2VF('?)(1/") 



hm , , , , 

n^oo 2W(<i){2/n 

1 1 



l/n2 



2T4^(9)'(0) 



— iy(«)"(o) < 
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and limn^aoiv{a+l/n) — v{a—l/n)) = by the continuity of the value function. 
Hence increment (1451) also tends to as n — > cxi. Taking limit as n — >■ c» in 
(|lll)-(|ini) completes the proof of (gH). □ 
All results of this section could be summarized in the following theorem. 

Theorem 5 The value function v{x) corresponding to the strategy iTa is given 
in lEBi) - 



6 Examples 

6.1 Classic risk process ([T]) with exponential jumps 

Assume that Xt is Cramer-Lundberg risk process ([T]) with exponential claims 
F{dz) — £^e^^^ dz and intensity of their arrival A. Then under P*^'?), where 



$(<z) 



q + X-Cc+ ^((7 + A-Cc)2 +Acq^ 



2c 



process Xt is again Cramer-Lundberg risk process ([T]) with exponential claims 
with parameter S,q = S,+^{q) and intensity of Poisson arrival equal to Xq = A^/^^ 
(see 0). 

We start form the strategy 7r° , where we pay dividend until the Parisian ruin 
time. The value function for the barrier strategy is given in (|28l) with 



y(9)(x) =e*(«)" ( 1 



XqD 



C^q-Xqil-D) 



(47) 



where 



D = 1 




<K+^i.)H-^I^{2t,/^)dt 



and Ii (x) is modified Bessel function of the first kind (see [9] and [10] for details). 
Then solving (|3ip we derive the optimal barrier: 



- A, 



log 



XqD 



C^q~Xq{l-D] 



1 - 



C^{q) 



From Corollary [2] it follows that barrier strategy with barrier a* is optimal 
strategy. 

In the second scenario tTq we have delay between decision to pay and its 
implementation. Then the value function is given in ()38|) with the scale function 
given by: 



where A± = q+^(^)-q-\q) with q'^{q) = ^(q) and q (q) being the smallest root 
of ip{0) = q: 

q + X~^C- y/{q + X-^c)^+Acq^ 



2c 



14 



Moreover, to identify v{a) note that, by ([5^ . ([T^ and lack of memory of expo- 
nential distribution, for a; > a we have, 

v{x) - E,_,e-«'^ [Xd, >d]+ v(a)e-9''P,_a (to" > d) 



+v{a)Ex-a 



l^(«)(c 



^e-^y dy 



y^x-a{T^ > d,Xde dy) 
+«(a)e-*'^P,_a {t- > d) 

Ma) I e-^^F.Mro € dt) J,, ^J^ e'^^ dy. (48) 



Further, 



nd nOG 

P,(to- >d,Xde dy) = P,{Xd e dy)- / P,(to- e dt) / P-^„(Xd_t e dy)^e-^ 

JQ Jo 

with P^X, e dy) = e-^«5cs(dy) + Er=i T^"^' (^^(c^-y)''"'^-^^^-^) dy. 
Finahy, by [4, Prop. IV.1.3, p. 101]: 

..(r-<t).Ae-(-VC).-ir^lif»d^, 



where 



/A 



A 



/i(6') = -exp<^ 2J-icos6i- (l + A/^)t + a;(J-cos6i- 1) ^ , 



'A 



A 



f2{9) = cos I — sin6' j — cos I — sin0 + 20 j , 



/3(0) 



1 + --2. 



i V c 



■ COS0. 



6.2 Brownian motion with drift - small claims 

Assume that 

Xt = aBt + ct, 

where cr, c > and Bt is a standard Brownian motion. Then under P**^?)^ where 
$(g) = with + 2q(T2^ have by © that = aBt + c^i. 

Considering strategy tt" with Parisian delay at ruin, the value function is 
given in ([25)1 with 
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where 



^'(a;) = 2^/ttxN{\/2x) - ^fnx + e" 



and is a cumulative distribution function for the standard normal distri- 

bution (for details see [5] and [H]). Hence by ([5T|) the optimal barrier is given 
by: 



= loe 

2c„ ^ 



' cr V 2 



and by Corollary [T] it is the optimal strategy. 

To get value function v{x) given in (j38]) for the case tTq with Parisian delay at 
the moment of payment of dividends, we have to identify first the scale function: 



1 



where 5 = a \J c? ^ Iqa'^ and a; = c/cr^. Then, similarly like in the previous 
example, 



v(x) 



-qd 



+v{a) 



yPx-a{To >d,Xd£ dy) 
e-«'^P,_a (r- >d)+ E,_, \e-^^o- , < dj 1 (50) 



and 



P,{to >d,Xde dy) = P,{Xd edy)- [ F,{t^ e dt)PiXd-t e dy) 

Jo 

with Xs having N(cs, s) distribution and Tq~ having inverse Gaussian distribu- 
tion: 



\{to e dt) 
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Dividend problem with Parisian delay for a 
spectrally negative Levy risk process 

Irmina Czarna* Zbigniew Palmowski^ 
October 19, 2011 



Abstract. In this paper we consider dividend problem for an insurance 
company whose risk evolves as a spectrally negative Levy process (in the 
absence of dividend payments) when Parisian delay is applied. The ob- 
jective function is given by the cumulative discounted dividends received 
until the moment of ruin when so-called barrier strategy is applied. Ad- 
ditionally we will consider two possibilities of delay. In the first scenario 
ruin happens when the surplus process stays below zero longer than fixed 
amount of time > 0. In the second case there is a time lag d between 
decision of paying dividends and its implementation. 
Keywords: Levy process, ruin probability, Parisian ruin, risk process, 
dividends. 



MSG 2000: 60J99, 93E20, 60G51. 

1 Introduction 

In risk theory we usually consider classical Cramcr-Lundberg risk process: 



where a: > denotes an initial reserve and Ui, {i — 1, 2, ...) are i.i.d distributed 
claims with the distribution function F. The arrival process is a homogeneous 
Poisson process Nt with intensity A. The premium income is modeled by a 
constant premium density c and the net profit condition is then XEUi/c < 1. 
Lately there has been considered more general setting of a spectrally negative 
Levy process. That is, X = {Xt}t>o is a process with stationary and indepen- 
dent increments having nonpositive jumps. We will assume that process starts 
from Xq = X and later we will use convention P(-|Xo = a:;) = and Pq = P. 

Such process takes into account not only large claims compensated by a steady 

'Department of Mathematics, University of Wroclaw, pi. Grunwaldzki 2/4, 50-384 
Wroclaw, Poland, e-mail: czarna@math.uni.wroc.pl 

tDepartment of Mathematics, University of Wroclaw, pi. Grunwaldzki 2/4, 50-384 
Wroclaw, Poland, e-mail: zbigniew.palmowski@gmail.com 




(1) 



income at rate c > 0, but also small perturbations coming from the Gaussian 
component and additionally (when v{—oo,0) — oo for the Levy measure v of 
X) compensated countable infinite number of small claims arriving over each 
finite time horizon. Working under this class of models, it became apparent 
that, despite of the diversity of possible probabilistic behaviors it allows to ex- 
press all results in a unifying manner via the g-harmonic scale function M^('J) (a;) 
defined via its Laplace transform. This paper further illustrates this aspect, by 
unveiling the way the scale functions intervenes in a quite complicated control 
problem. 

The classic research of the Scandinavian school had focused on determining 
the "ruin probability" of the process ^ ever becoming negative, under the 
assumption that X has positive profits. Since however in this case the surplus 
has the unrealistic property that it converges to infinity with probability one, 
De Finetti [14) introduced the dividend barrier model, in which all surpluses 
above a given level are transferred (subject to a discount rate) to a beneficiary, 
and raised the question of optimizing this barrier. Formally, we consider the 
risk process controlled by the dividend policy tt given by 

U^^Xt-L^, (2) 

where Xq = a; > is an initial reserves and is an increasing, adapted and 
left-continuous process representing the cumulative dividends paid out by the 
company up till time t. The optimization objective function is given by the 
average cumulative discounted dividends received until the moment of ruin: 

v^x)^E^ f e-«*dL^ (3) 
Jq 

where is a ruin time that we specify later depending on the considered 
scenario and g is a discounting rate. 

The objective of beneficiaries of an insurance company is to maximize v'^{x) 
over all admissible strategies tt: 

v*{x) — snpv^{x), (4) 
Tren 

where H is a set of all admissible strategies, that is strategies vr = {L^,t > 0} 
such that - < C/f_. 

For classical risk process ([Ij an intricate "bands strategy" solution was dis- 
covered by Gerber [TSJIH], as well as the fact that for exponential claims, this 
reduces to a simple barrier strategy: "pay all you can above a fixed constant 
barrier a" . 

There has been a great deal of work on De Finetti's objective, usually con- 
cerning barrier strategies. Gerber and Shiu [17] and Jeanblanc and Shiryaev [21] 
consider the optimal dividend problem in a Brownian setting. Irback [50] and 
Zhou [33] study the constant barrier under the Cramer-Lundberg model ([1]). 
Hallin [19| formulated time dependent integro-differential equations describing 
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the payoff associated to a 2n bands policy. The optimahty of the "bands strat- 
egy" was recently established by Albrecher and Thonliauser [1] in the presence 
of fixed interest rates as well. For related work considering both excess-of-loss 
reinsurance and dividend distribution policies (e.g. in a diffusion setting), see 
Asmussen et al. [3] and |29| and references included in this papers, and for work 
including also a utility function, see Grandits et al. [TH]. 

For Levy risk process considered in this paper without Parisian delay Avram 
et al. [B], LoeSen and Loeffen and Renaud [17] found sufficient conditions 
for which barrier strategy is optimal. In fact Avram et al. [7] prove that in this 
case the bands strategy is optimal. 

In this paper we want to analyze the dividend problem when so-called 
Parisian delay either at moment of dividends payments or at the ruin time 
is applied. The name for this delay comes from Parisian option that price are 
activated or canceled depending on the type of option if the underlying asset 
stays above or below the barrier long enough in a row (see [5] and [13] )■ 

For the classical risk process ([T]) Dassios and Wu [H] consider a Parisian 
type delay between a decision to pay a dividend and its implementation. The 
decision to pay is taken when the surplus reaches the fixed barrier a but it is 
implemented only when the surplus stays above barrier longer than fixed d > 0. 
The dividend is paid at the end of this period. This strategy we will denote 
by TTa,d- Similar problem for a spectrally negative Levy process of bounded 
variation was analyzed in [25_. In this paper we generalize this result into 
general spectrally negative Levy risk process. In this case ruin time is given 
by: a"""-" = inf{t > : Ut"'" < 0}. Since the ruin time is classical wc know 
that band strategy is optimal and we also know the necessary conditions under 
which an optimal strategy is the barrier strategy. We still believe that this 
new Parisian strategy (although not optimal within all strategies) could be very 
useful for the insurance companies giving possibility of natural delay between 
decision and its implementation. 

In this paper we also consider Parisian delay at the ruin. We denote this 
strategy by n'^'^ . That is ruin occurs if process stays below zero for longer 

period than a fixed C > 0. Formally, we define Parisian time of ruin by: 

cr''° '^ = mi{t >0:t- sup{s < t : Uf'^ > 0} > C, Uf'^ < 0}. (5) 

Wc first analyze the strategy tt"'^ for which the dividends are paid according to 
classical barrier dividend strategy transferring all surpluses above a given level 
a to dividends. We also prove the verification theorem for this type of ruin. In 
particular we find sufficient condition for the barrier strategy to be optimal. 

In fact combination of both scenarios is also available. To simplify analysis 
we decided to skip this possibility. 

We believe that giving possibility of Parisian delay could describe better 
many situations of insurance company. For example it can be checked if indeed 
company's reserves increase and we can pay dividends (in the first scenario) or 
it can be given possibility for the insurance company to get solvency (in the 
second scenario). 
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The paper is organized as follows. In Section [3] we introduce basic notions 
and notations. In Section [3] we find discounted cumulative dividends payments 
until Parisian ruin time. In Section we prove the verification theorem and 
find necessary conditions for the barrier strategy to be optimal. In Section [5] 
we analyze the case when there is a time lag between decision to pay dividends 
and its implementation. Section [S] is devoted to examples. 



2 Preliminaries 

We first review some fluctuation theory of spectrally negative Levy processes 
and refer the reader for more background to Kyprianou [53], Sato [32] and 
Bertoin [8] and references therein. 

In this paper we consider a spectrally negative Levy process X — {Xt\t>o, 
that is a Levy process with the Levy measure v satisfying v{0, oo) — (for 
simplicity we exclude the case of a compound Poisson process with negative 
jumps) . Since jumps of a spectrally negative Levy process X are all non-positive, 
moment generating function E[e^'^*] exists for all > and is given by E[e^'^*] = 
Qt^(s) fQj. gome function ipifi) which is strictly convex with the property that 
lim^-^oo V'(^) = +00. Moreover, ij; is strictly increasing on [$(0), oo), where $(0) 
is the largest root of ^(6') = 0. We shall denote the right-inverse function of -0 by 
$ : [0, oo) — >■ [$(0), oo). We will consider also the dual process Xt = —Xt which 
is a spectrally positive Levy process with the Levy measure V (0, y) = v {—y, 0). 
Characteristics of X will be indicated by using a hat over the existing notation 
for characteristics of X. 

For any 9 for which ■ip{d) ~ \ogE[exp 9Xi] is finite we denote by an 
exponential tilting of measure P with Radon-Nikodym derivative with respect 
to P given by 

^ =exp{eX,-mt), (6) 

t 

where Ft is a right-continuous natural filtration of X. Under the measure P^ the 
process X is still a spectrally negative Levy process with characteristic function 
ifje given by 

^e(s) =V(s + 0)-V'(^). (7) 

Throughout the paper we assume that the following (regularity) condition 
is satisfied: 

.0 

cr > or / xvidx) = oo or v{dx) « dx, (8) 
where a a Gaussian coefficient of X. 
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2.1 Scale functions 



For p > 0, there exists a function W^^^ : [0,oo) — >■ [0,oo), called the p-scale 
function, that is continuous and increasing with Laplace transform 



-Ox 



e > <s>{p). 



(9) 



The domain of (p' is extended to the entire real axis by setting W'^p^ (y) = for 
2/ < 0. We denote W'-^^z) = W{z). For later use we mention some properties of 
the function W^^^ that have been obtained in literature. On (0, oo) the function 
z H-)- W^P\z) is right- and left-differentiable and under the condition (|5]), it holds 
that z I-)- W^P^y) is continuously differentiable for y > 0. Moreover, if cr > 
it holds that W^p'> e C°°(0,oo) with W^p'>'{0) = if X has unbounded 

variation with cr = 0, it holds that W'^p^'{0) = oo (see [30] Lemma 4]). 

The function plays a key role in the solution of the two-sided exit 

problem as shown by the following classical identity. Letting t+ , t~ be the 
entrance times of X into [a, oo) and (~oo, —a) respectively, 



T+ = inf{t >0:Xt>a}, 
it holds for z £ [0, a] that 

E 



6 tT^O > 



M{t >0:Xt< -a} 



= W^P\z)/W'^P\a). 



(10) 



(11) 



Closely related to W'^p^ is function Z^p^ given by 

Z(f)(z) = l + gW^^'^z), 

where w''^\z) — W'^P\y)dy is the anti-derivative of W^pK Moreover, the 
scale functions appear also in so-called two-sided downward exit problem: 

M/(p)(z) 



^ ' n 



Z^P\z) - Z'^P\a)- 



(12) 



M^(p)(a)' 

and in one-sided downward exit problem that for any /3 with ■0(/3) < oo> V ^ 
t/'(/3) V and a; > gives: 



where W^"'' and Z^"' are scale functions with respect to the measure P^, u = 
p — ?A(/3) and ^/^(u) is understood in the limiting sense if u = 0. In fact for 
each z G R, W^'p^ (z) is analytically extendable, as a function in p, to the whole 
complex plane; and hence the same is true of Z^p\z). In which case arguing 
again by analytic extension one may weaken the requirement that p > tpil^) V 
to simply p > 0. 

The 'tilted' scale functions can be linked to non-tilted scale functions via the 
(z) = W^P'){z) from O Remark 4]. This relation implies that 



-PT-+I3X _ 

e 



An) 



Tg < oo 



.0^ 



z'-\z) 



$(w) 



(13) 



relation e^^w'"^^ 



Z^^\z) = 1 + u 



-PyW^P\y)dy. 
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2.2 Parisian ruin 



One of most important characteristics in risk theory is a ruin probability defined 
by Vxi^Q < oo) for Tq" = mi{t > . Xf < 0}. Czarna and Palmowski [5] 
extended this notion to so-called Parisian ruin probability, that occurs if the 
process X stays below zero for period longer than a fixed C > (see also [TUl H] 
for the result concerning classical risk process). Let 

= inf {i >0:t~ sup{s < t : Xs > 0} > (, Xt > 0} 

and Parisian ruin probability we define as: 

P(r'^ < oo\Xn = x)= P^ir'^ < oo). 

The following result summarize Theorem 1]. 

Theorem 1 For any x > Parisian ruin probability equals: 

P^(r'- =00)= E[Xi]-^ r'^ \ ' 

Jo ^IP(^C ^ dz) 

3 Parisian delay at ruin 

In section we will consider Parisian ruin time ([5|) and dividends paid according 
to barrier strategy that correspond to reducing the risk process to the 

level a if x > a, by paying out the amount x — a, and subsequently paying out 
the minimal amount of dividends to keep the risk process below the level a. It is 
well known (see [B]) that for < x < a the corresponding controlled risk process 
U"^ ^ under P^; is equal in law to the process {a — Yt : t >0} under for 

Yt^aWXt-Xt (14) 

being Levy process X reflected at its past supremum: 

Xt = sup Xs, 

0<s<t 

where we use notations y V = ma.x{y, 0}. In this case for all a; > 0, 




and ' = aV Xt — a. 
Note that for x < a, 



v''''^{x)=E, e-'i^^,T+ <t'^ w"'«(a) (15) 
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and 

v'^''^ {x) ^ X - a + v"''^ {a) for x > a. (16) 

Assume that X oo a.s. Then by Markov property and fact that X jumps 
only downwards we derive: 

P^(r« = oo) = P,(t+ < ri)¥a{T'^ = cx)). (17) 

Hence 

Pa(r? ^ oo) 

Using change of measure ^ with 9 ~ <&(?), Optional Stopping Theorem and fact 
that on P*^'?) process X tends to infinity a.s. (since V'$(q)(0+) = V'''(*^'('7)+) > 
0), we have for x < a, 

V^i\x) 



where by Theorem [U 

F(9)(x) = e*(«)"P*(«)(T^ = oo) = E*(«)[Xi] 



- Vi^){a)' ^^^^ 
e~Hg)-W^g)(^x + z)zP*(9)(X^ e dz) 



/p°°zP*(9)(X(^ € dz) 



(19) 



It follows from Theorem [T] that under the condition ([SJ function F^'^ (y) 
(similarly like W^'^^{y)) is continuously differentiable for j/ g R. 
Moreover, for n € N, by (HI 



' a+l/n 



f"'«(a) + - 
n 



and 



1 

— E„ 
n 



1 - E, 

nq 



' a+l/n 

(^w'^'^(a) + 

' a+l/n y 



+ T-+ < T-C 



since LJ^ ' = Xt — a under P^ can increase only by e up to time t^_^_i/j^- Last 
increment in above equation is o(l/n) since by strictly positive drift X is regular 
for (0, oo). Hence, 
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and then 



^ ' V('i)'{a) 



Thus from p^ . (IT5|) and it foUows that is continuously differentiable 
for ah a; e M and 

(20) 

In particular, 

= 1- (21) 

Hence we get the following theorem. 

Theorem 2 The value function corresponding to the barrier strategy 7r"'^ is 
given by i20\). The optimal barrier a* is given by: 

a* = inf{a > : V'^''^'{a) < V'^'^^'iy) for all y > 0}. (22) 

In particular, if V^'^'^ € C^(M) and there exists unique solution of equation: 

V^^^"{a*) ^ 0, (23) 

then a* is optimal barrier. 

Remark 1 Note that € C^iR) if W^i'^ G C^iR). This is the case if e.g. the 
Gaussian component is present. 



4 Verification Theorem 

To prove the optimality of a particular strategy tt across all admissible strategies 
n for the dividend problem where the ruin time cr^ is given by the Parisian 
ruin ([5]), we are led, by standard Markovian arguments, to consider the following 
variational inequalities: 

Tf{x)-qfix) < 0, if x€M, (24) 
fix) > 1, if X e R, (25) 

for functions / : R — > R in the domain of the extended generator F of the process 
X which acts on C^(0, oo) functions / as 

r/(x) = y/"(a;) +Ponx) + J" [fix + y)- fix) + /'(x)yl{|,|<ij] i^idy), 

(26) 
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where v is the Levy measure of X and cr^ denotes the Gaussian coefficient and 
= c — l^_^yv{dy) if the jump-part has bounded variation; see Sato [3H 
Ch. 6, Thm. 31.5]. In particular, if E\X\ < oo and X has unbounded variation 
[bounded variation], a function / that is [C^] on [0, oo) and that is uhiniatcly 
hnear hes in the domain of the extended generator. 

Theorem 3 Let C £ (0,oo] and suppose f is continuous and piecewise on 
(— oo,C) if X has hounded variation and that f is and piecewise on 
(— oo,C) if X has unbounded variation. Suppose that f satisfies \2^^ and i25\) . 
Then f > sup^^jj^^ v'^ for defined in with Parisian ruin time 0), where 
n<c is a set of all bounded strategies by C. In particular, if C = oo, then 

Proof We will follow classical arguments. Let tt G Il<c be any admissible 
policy and denote hy L = and U = the corresponding cumulative div- 
idend process and risk process, respectively. By Sato [32l Ch. 6, Thm. 31.5] 
function g{t,x, z) = e~'^^ f {x) \ ^^kq is in a domain of extended generator of the 
three-dimensional Markov process (i, , <,y), with = i — sup{s < t : ?7t > 0}. 
Note that finite number of discontinuities in / and hence also single discontinu- 
ity in \^z<Q} are allowed here. Hence we are also allowed to apply Ito's lemma 
(e.g. [311 Thm. 32]) if X is of unbounded variation and the change of variable 
formula (e.g. [31] Thm. 31]) if X is of bounded variation: 

e-«*/(t/*)l{.f<c} - /(C^o) = J}{t)- f e-'^^f'(U,-)dL'; 

Jo 

+ f e"^'{Tf~qf){Us-)ds + Mt, (27) 

JQ 

where Aft is a local martingale with Mq — 0, L'^ is the pathwise continuous 
parts of L and for a function g the process Jg is given by 

Ut) = E e''' + Bs) - g{A,)] l{B.^o}, (28) 

s<t 

where Ag — Ug^ + AX^ with Axg = Xg — Xg- and Bg = —ALg denotes 
the jump of —L at time s. Let r„ be a localizing sequence of M. Applying 
Optional Stopping Theorem to the stopping times T^. = A cr^ and using 
Fatou's Theorem we derive: 

fix) > E,e-''^"/(C/T;)l{,^, <c} - Jf{K)+E, f " e~'^'f'iUg-)dL'g 

Jo 

-E, r" e"'^rf-qf){Ug-)ds. 
Jo 
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Invoking the variational inequalities /'(x) > 1 (hence f{As + Bg) — f{As) < 
-ALs if As > 0) and Tf{x) - qf{x) < we have: 



fix) > E,e-«^"/(C/Tjl{,c.,<c}+E, 











Jo 








+ E, 


r e-'^dL,;(T- 




Jo 



Letting n — oo in conjunction with the monotone convergence theorem and 



using fact that 1|^!7^<^| — complete the proof. 



□ 

Using verification theorem we find necessary conditions under which the 
optimal strategy takes the form of a barrier strategy. 

Theorem 4 Assume that a > Q or that X has bounded variation or, otherwise, 
suppose that v'^ •'' G C^(0,oo). If q > 0, then a* < oo and the following hold 
true: 

(i) it"' is the optimal strategy in the set Il<a' of all hounded strategies by 
a* and v"" = sup^gn^^, . 

(ii) If {Vv"" '^ - qv"' 

)i^) ^ for X > a* , the value function and optimal 
strategy of is given by = v" where the ruin time cr^ is given by the 
Parisian moment of ruin 

The proof of Theorem |4] is based on the verification Theorem |3] and the 
following lemma. 



Lemma 1 (i) We have a* < oo. 

(ii) It holds that (Lu"*'? - gw"*'«)(a;) = for x < a*. 
(Hi) For X < a* 

> 1. 

Proof Point (i) follows from the fact that V"''?''(j/) is continuous and in- 
creasing starting at some point. Indeed, note that by [22 we have V'-'^'^iy) = 
e'f{q)ypyir< = oo) > e'^^'^^yPyir^ = oo) ^ W^^^'^^lv) ^nd tends to 

00 as ?; — >■ oo. The proof of (ii) follows from ()15p and the martingale property 
of 



Xt 



I ' a* ^ ' 



E, 



I 'a* ^ ' 



where x < a*. Point (iii) is a consequence of (|2ip and definition of a* given in 
(EH). □ 
Moreover, we can give other necessary condition for the barrier strategy to 
be optimal. 

Corollary 1 Suppose that 

y('?)'(a) < for alia* <a<b. 

Then the barrier strategy at a* is an optimal strategy. 
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Proof Using Theorems [5] and S] the proof is the same as the proof of |26[ Theorem 
2]. □ 



Corollary 2 Suppose that, for x > 0, v'{x) is monotone decreasing, then tt" 
is an optimal strategy of 

Proof By the proof is similar like the proof of Theorem 3] . In fact it 
suffices to prove that V$(q) has completely monotone derivative. This fact follows 
from Theorem [U since ^Ff-^^T^ = oo) and ^P*^''^(T|7 < oo, -X^- e dz) are 
completely monotone. Indeed, it is known that 

^P*(«)(ro- =^*(,)(0,0)^4(^)(0,a:), (29) 

where J7*(g) is a renewal function of the descending ladder height process Ht 
under P*(9) and (a, /?) is a Laplace exponent of bivariate descending ladder 
height process (L^^^Ht) under P**?) with k<e,(,)(0,0) = i}'{<^{q)) > 0. From the 
proof of [221 Theorem 3] it follows that C^4(g)(^' ^) hence also ^P*^'^(rg' — oo) 
is completely monotone. Moreover, by [24l (7.15), p. 195] we have 

= K$(9)(0,0) / U^^g){x-dy) e-'^^^^^'+^^V{dz + v) dv, 
Jo Jo 

and hence ^P*^'^'(to" < oo, —X^- e dz) is also completely monotone. □ 



5 Parisian delay at the moment of dividend pay- 
ments 

In this section we analyze the case when we pay dividends only when surplus 
process stay above barrier a longer than a time lag d > 0. The dividends 
are paid at the end of that period and they are paid until regular ruin time 
= M{t > : Ut''-'' < 0}. Then by ^ iov x € [0, a], 



'"aA"-) = wr.w N ^a,d(a); (30) 



1^(9) (a) 



and by Markov property for x > a, 



(0,a) 



« , -X- £ dy, To 



11 



where E,. 



I 'a ^ 'O 



Double Laplace transform of 
T3|) equals 

roc r 



is given in ([TU|) . 



"'[0,00) 



for z > by 



ds 



-ia+q)T-+l3X _ _ 

e ^0 , Tq < 00 



(32) 



where Uq ^ a + q — Moreover, 



X -e "^0 , Tn < 00 



iz-^i/o(/3,^)|,=o. 



(33) 



Further, the value Wa,d(a) is determined by pip if X has no Gaussian com- 
ponent (ct = 0) or by the smooth paste condition: 



(34) 



otherwise. 



Lemma 5.1 If a > then holds. 
Proof For n G N, 



+Wa,<i(a + l/^^)Ea 



' 'a-l/n ^ 'a+1/7^ 



" nA-\ IX n 



a+l/?i a— 1/n 



(35) 



where the last term is bounded above by -^PiT^^^ '^)- Moreover, by pH)) 



14^(9) (1/n) 
W(i'){2/n) 



and by ([12 

E, 



' a— 1/n a+l/n. 



Z('')(l/n)-Z(«)(2/n) 



M^(9)(2/n) 
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Multiplying both sides of ([551) by two, subtracting VaAia—^ /n)+Va,d{a), dividing 
by 1/n gives: 



fa,d(a) - fa,rf(a - 1/"-) 
1/n 



VaA^ + 1/") ~ •^a,d(g) 

1/n 

VaA{a + 1/n) - Wa,d(a - 1/n) 



l/n 



2M^('?)(l/n) 
H^(9)(2/n) ^ 



(36) 



(37) 



+Wa,d(a - l/n)2q 



l/n 



where we use (|lip . Now, since 

W^'HO) = and 1^(9)' (0) = 4, we have: 
VK(«)(l/n) M^(9)'(0) _ 1 



lim 



W^i^{2/n) 2W(9)'(0) 



Hence increment ^ converges to Wa,d(a)2(?(VF('?nO)-^Vl/(«)(0)) = as n — oo. 
Moreover, 



lim n 



2W^(«)(l/n) 
W('i){2/n) 



1 2/n H/(9)(2/n) - 2T/F('?ni/") 
lim 



n^oo 2V7(9)(2/n) 



l/n2 



and lim„^oo(z;a_d(a + 1/n) — Va,d{o- — 1/n)) — hy the continuity of the value 
function. Hence increment (1571) also tends to as n — > oo. Taking limit as 
n oo in (l36|)-(l38l) completes the proof of □ 
All results of this section could be summarized in the following theorem. 

Theorem 5 The value function Va,d{x) corresponding to the strategy TTa.d is 
given m ^ - 



(38) 



6 Examples 

6.1 Classic risk process with exponential jumps 

Assume that Xt is Cramer-Lundberg risk process ([1]) with exponential claims 
F{dz) = S,e~^^ dz and intensity of their arrival A. Then under P*^'?), where 

. g + A - + v/(g + A - $c)2 + Acqs, 

m = Yc ' 

process Xt is again Cramer-Lundberg risk process ([1} with exponential claims 
with parameter ~ C+"^(9) and intensity of Poisson arrival equal to \ = ^i/iq 
(see 0). 
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6.1.1 Dividend strategy tt" '' 



We start form the strategy 7r°''», where we pay dividend until the Parisian ruin 
time. The value function for the barrier strategy is given in (1201) with 

W..W...<..-(l-e-(^)-(-^_M^)). (33, 

where 

Jo y 

and Ii (x) is modified Bessel function of the first kind (see [5] and [TU] for details). 
Then solving we derive the optimal barrier: 

c^,^\,{l-D,)) \ c^{q) ) _ ■ 

From Corollary [5] it follows that barrier strategy with barrier a* is optimal 
strategy. 

We present now the numerical analysis. We study Parisian dividends for 
different Parisian delays (Table 1): ^ = 0.1, 0.3, 0.7, 2 and for different initial 
capitals (Table 2): a; = 2, 5, 10, 50. Moreover in Table 3 we compare the initial 
capitals when dividend payments are made till classical and Parisian ruin time. 
In the first case we denote the value function by v (see [6] for details concerning 
identifying it). We choose the following parameters: ^ = 2, A = 2, c = 2.5 (see 
also Albrecher and Thonhauser [1]). 



loa 



c 


0.1 


0.3 


0.7 


2 


a* 


3.54 


3.09 


2.40 


0.84 



Table 1: The optimal barrier for different Parisian delays. 



X 


2 


5 


10 


50 


v{x) 


12.57 


15.71 


20.71 


60.71 




13.38 


16.40 


21.40 


61.40 



Table 2: The value function for classical and Parisian dividend problem for 
different initial capitals, when Parisian delay equals ^ = 0.3. 



X 


2 


5 


10 


50 


Xl 


2.69 


5.69 


10.69 


50.69 


v{xi) = v°-' '''{x) 


13.38 


16.40 


21.40 


61.40 



Table 3: The same value function for classical and Parisian dividend problem 
for different initial capitals, when Parisian delay equals Q — 0.3. 
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6.1.2 Dividend strategy iTa.d 

In the second scenario iTa^d we have delay between decision to pay and its im- 
plementation. Then the value function is given in pO[) with the scale function 
1^(9) given by: 



where A± = q+^(^^-q-\q'f with q'^{q) = ^{q) and q (q) being the smallest root 
of - q- 



_ g + A-gc- v/(g + A-^c)2+4^ 

('^^ = 2-c • 

Moreover, to identify Va^dia) note that, by ([5T|) . and lack of memory of 
exponential distribution, for a; > a we have. 



VaA^) = E,_„e-«'^ [Xd, > d] + VaAa)e-'"'TP^-a (tq- > d) 
-|-WQ,d(a)E^_Q 



M/(?)(a) ^ ^ 



Jo 

+«a,d(a) e-^*P,_,(ro- e dt) £ ^^(1) (~)^^ ^^'^^ dy. (40) 



Further, 



/•d POO 

.(tq- > d, e dy) - P,{Xd e dy)- / P3(to- G d<) / P-i„(Xd_t G dy)S, 

Jo Jo 



e «™ 



withP,(X, G dy) = e-^^5e,(dy) + Er=i ^e^^^ (^(cs-y)^-ie-«(— dy. 
Finally by [4l, Prop. IV.1.3, p. 101]: 



where 



fi(.0) = ^exp|2y^icos0-(l + A/e)t + a;(y^cos0-l)|, 
f2{9) = COS ^^J^sinO^ — cos + , 

fsiO) = l + ^-2J^cos0. 
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6.2 Brownian motion with drift - small claims 



Assume that 



ct. 



where a,c > and is a standard Brownian motion. Then under P*^'?) where 



with Cq — ^/cP + 2qa^, we have by ([7]) that Xt — aBt 



Cqt. 



aX 



6.2.1 Dividend strategy tt 

Considering strategy tt"^'^ with Parisian delay at ruin, the value function is given 
in (HOl with 



^<S>(q)x 



(41) 



where 



*(a;) = 2y/7TxJ\f{V2x) - V^i^x + 



and A/'(.) is a cumulative distribution function for the standard normal distri- 
bution (for details see [5] and [H]). Hence by the optimal barrier is given 
by: 



£2. 



1 - 



2c„ 



and by Corollary [T] it is the optimal strategy. 

Similarly like for the Cramer-Lundberg process we present here the numerical 
analysis. We study Parisian dividends for different Parisian delays (Table 4): 
C = 0.1, 0.3, 0.7, 2 and for different initial capitals (Table 5): x = 2, 5, 10, 
50. In Table 6 we compare initial capitals for classical and Parisian dividend 
problems. We take c = 2.5 and cr = 2. 



c 


0.1 


0.3 


0.7 


2 


a* 


4.48 


3.89 


3.12 


1.17 



Table 4: The optimal barrier for different Parisian delays. 



X 


2 


5 


10 


50 


v{x) 


20.49 


24.72 


29.72 


69.72 


v'^'^'^ix) 


23.00 


26.11 


31.11 


71.11 



Table 5: The value function for classical and Parisian dividend problem, when 
Parisian delay equals ( = 0.3. 
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X 


2 


5 


10 


50 


Xi 


3.40 


6.39 


11.39 


51.39 


v{xi) = i'"'''»{a:) 


23.00 


26.11 


31.11 


71.11 



Table 6: The same value function for classical and Parisian dividend problem, 
when Parisian delay equals C = 0.3. 



6.2.2 Dividend strategy TTa.d 

To get value function v{x) given in (pO)) for the case iTa.d with Parisian delay at 
the moment of payment of dividends, we have to identify first the scale function: 



where S = a ^-\/c^ + 2qa^ and ui = c/a^. Then, similarly like in the previous 
example, 



-qd 



yPx-a{TQ >d,Xd£ dy) 



< d 



(42) 



and 



P,(ro- >d,Xde dy) - P^iXd edy)- / P,{t^ G dt)PiXd-t G dy) 

Jo 

with Xs having N(cs, s) distribution and Tq having inverse Gaussian distribu- 
tion: 

Px{tq e dt) = 
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